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Abstract 

We prove uniqueness and Monge solution results for multi-marginal 
optimal transportation problems with a certain class of surplus functions; 
this class arises naturally in multi-agent matching problems in economics. 
This result generalizes a seminal result of Gangbo and Swiech on multi- 
marginal problems. Of particular interest, we show that this also yields a 
partial generalization of the Gangbo-Swi§ch result to manifolds; alterna- 
tively, we we can think of this as a partial extension of McCann's theorem 
for quadratic costs on manifolds to the multi-marginal setting. 

We also show that the class of surplus functions considered here neither 
contains, nor is contained in, the class of surpluses studied in ^28J, another 
generalization of Gangbo and Switch's result. 

1 Introduction 

Two marginal optimal transportation is the general problem of coupling two 
distributions of mass together as efficiently as possible, relative to a given surplus 
function. This is an exciting and very active area of research, with a wide variety 
of applications. For a full literature review, see the book of Villani |32j . 

Recently, optimal transport problems with several marginals have begun 
to attract more attention, due to emerging applications in economics [S] [11] . 
mathematical finance [IS] [12] [1] , condensed matter physics [H] [5] , barycenters 
and image processing [I] [3] and the characterization of m-cyclically monotone 
vector fields [14] . This is the general problem of aligning several mass dis- 
tributions with maximal efficiency. Stated precisely, given Borel probability 
measures fn on bounded, open sets Mi C M", for i = l,2...m and a surplus 
function b : Mi x M2 x ... x Mm — > K, the problem is to maximize 
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/ b{xi,X2,...,Xm)dj (MK) 

JAfixAf2X...xM,„ 

among measures 7 on Mi x M2 x ... x Mm which project to the ^i. Under 
reasonable conditions, a maximizer 7 exists. We say that 7 is a Monge solution 
if it is concentrated on the graph of a function over the first variable. 

When m = 2, this is the classical, two marginal Monge-Kantorovich problem. 
In that setting, assuming the first marginal is absolutely continuous with respect 
to Lebesgue measure and the surplus function b satisfies a fairly weak condition, 
called the twist condition, it is by now well known that the solution 7 induces 
a Monge solution and is unique [H] [IS] [7] . 

When m > 2, however, the picture is much muddier. In contrast to the two 
marginal case, the structure of solutions to these problems is not well under- 
stood, except in certain special cases. Uniqueness and Monge solution results 
have been proven by Gangbo and Swi§ch [17] for the quadratic surplus, 

m 

h[xi,X2, —,Xrn) = - ^ l^;* - Xjf; (1) 

see also [24] [20] [30] [31] . This was generalized by Heinich [18] to surplus functions 
of the form 

m 

b{xi,X2,.:,Xm) ^ h(^Xi), (2) 

1=1 

where h : M" M is strictly convex. When n — 1, Carlier was able to further 
generalize this result to the class of strictly 2-monotone surplus functions. In 
higher dimensions, a result of the present author establishes uniqueness and 
Monge solutions under strong second order conditions on the surplus [28], en- 
compassing the results in [17] and [18]. In contrast to the two marginal case, 
these results depend strongly on the form of the surplus function; a series of 
counterexamples in [10] [29] and [27] indicate that solutions may not be unique 
or induce Monge solutions, even when the surplus functions satisfy naive multi- 
marginal generalizations of the twist condition. In fact, the examples in |29j are 
supported on submanifolds S C Mi x M2 x ... x Mm of dimension n{m — 1). 

In this paper we study a multi-marginal optimal transportation problem for 
a particular class of surplus function. Given some other open subset Z C M" 
and functions fi : Mi x Z — > K, we will consider surplus functions of the form: 

m 

b{xi,X2,-.Xm) = supy^fi{xi,z) (3) 

Surplus functions of this form arise naturally in multi-agent matching, or 
hedonic pricing, problems in mathematical economics. These types of prob- 
lems were originally introduced by Ekeland 13 (when m = 2) and Carlier and 
Ekeland [9] (when to > 3); the formulation as an optimal transport problem 
with a surplus of the form ([3]) can be found in the paper of Chiappori, McCann 
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and Nesheim [TT]. Economically, each measure fii represents the distribution 
of a category of agents on the type space Mi (parameterizing, for example, the 
agents' skill sets). These distributions must be matched to make teams; each 
team consists of one agent from each category. The matching is done via signing 
contracts z G Z and the functions fi{xi, z) represents the preference of an agent 
of type Xi for contract z. The result of Chiappori, McCann and Nesheim tells us 
that finding an equilibrium in this market is equivalent to solving the optimal 
transport problem with surplus ; the optimal measure 7 tells us which agents 
come together to form teams in equilibrium. 

Our main goal in this paper is to resolve Monge solution and uniqueness 
questions for this class of surplus functions. Aside from its importance in eco- 
nomic applications, this class of surplus functions includes those studied by 
Gangbo and Swi§ch ([T]) and Heinich ([2]) (as we show in subsections 5.1 and 
5.2), and so our main theorem can be seen as a generalization of their results. 
Let us recall that the class of surplus functions treated in [55] also included ([T]) 
and (HI). As we show in the subsections 5.3 and 5.4, our result here is neither 
strictly stronger nor weaker than the result in |28| : that is, we exhibit explicit 
examples of surplus functions satisfying the hypotheses in |28| which are not of 
the form ([3]), as well as examples of the form ([3]) which do not satisfy the as- 
sumptions in [28 . It is interesting to note that our assumptions on the functions 
fi (conditions HI -H5 in the next section) look much more like typical optimal 
transportation type conditions than the strong, complicated looking conditions 
required in [28] . 

Unlike the results in [17] [l^ [8] and [28] , our argument here does not rely on 
a dual formulation of (|MKp . Instead, the proof of our main result (Theorem 
14. ip exploits two main ingredients. The first is a general, structural result on 
solutions to multi- marginal problems [53]. For surplus functions of the form 
this result surprisingly implies that the support of the solution n-rectifiable and 
spacelike for a certain semi-Riemannian metric. 

The second ingredient is the interplay between (jMKp and another variational 
problem, introduced by Carlier and Ekeland [5] (see (jMAMp in the next section). 
From an economic perspective, this second problem can be seen as another 
formulation of the multi-agent problem. In this formulation, roughly speaking, 
one focuses on the couplings between each category of agents and the contracts 
they sign (rather than on the coupling between the various categories of agents, 
as in the formulation of Chiappori, McCann and Nesheim). The problem is to 
find the probability measure v on the space of contracts which maximizes the 
sum of the two marginal transportation surpluses with each /z^; economically, 
this measure tells us the relative frequency of contracts that get executed in 
equilibrium and the optimal couplings between each /i^ and v tell us which 
types in each category sign which contracts. 

Roughly speaking, the rectifiability result from [53] allow us to apply tools 
from geometric measure theory (namely, the co-area formula) to the correspon- 
dence between teams and contracts (stated precisely, this is a map z : spt{'-f) — >■ 
Z, pushing 7 forward to ly). We then use the spacelike condition to show that 
different teams {xi,X2, ...,Xm) in the support of the optimizer always sign dif- 
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ferent contracts, at least locally, allowing us to deduce that v is absolutely 
continuous with respect to Lebesgue measure. From here, we can use standard, 
two marginal optimal transport techniques to establish that in equilibrium, un- 
der appropriate assumptions, almost all contracts z ^ Z are signed by at most 
one agent type Xi := Fi{z) S Mi (which easily yields the uniqueness of 7), and 
that the map Fi : Mi Z \s invertible almost everywhere. We then have that, 
for (xi, a;2, Xm) in the support of the optimizer, Xi = Fi(F^^(xi)), implying 
that 7 induces a Monge solution. 

As a consequence of our work here (in fact, as a direct ingredient in our 
argument) we prove a regularity result on the optimizer in the matching problem 
formulation of Carlier and Ekeland. In addition to its economic implications, 
this corollary may have independent interest. In the special case when each 
fi{xi,z) — — \xi — z\'^, the Carlier-Ekeland formulation is equivalent to finding 
the barycenter of the measures /x^, a problem introduced and solved by Agueh 
and Carlier [T]. Our regularity result generalizes a result by Agueh and Carlier, 
who use a correspondence between the multi-marginal problem with quadratic 
cost ([T]) and the problem of finding a barycenter of several measures on R"; 
in that paper, they used this relationship to establish an L°° estimate on the 
barycenter (assuming that at least one of the has an L°° density). At the 
same time, our approach can be seen as generalizing an argument of the present 
author in the two marginal case; in that work, the geometry of the optimal plan 
is used to establish regularity of v [55]. 

It is also natural to consider the barycenter problem in the Wasserstein 
space over a Riemannian manifold. Riemannian manifolds with non-positive 
sectional curvature fit naturally into the framework of this paper. We will show 
that our result also implies a regularity result on the barycenter of measures 
on simply connected, non-positive ly curved Riemannian manifolds, generalizing 
the regularity result of Agueh and Carlier on Euclidean space. 

Furthermore, in this setting, we show that our main result can be interpreted 
as a generalization of the theorem of Gangbo and Swi§ch to simply connected, 
negatively curved Riemannian manifolds, in the same spirit that McCann's po- 
lar factorization theorem on Riemannian manifolds generalizes Brenier's polar 
factorization theorem on Euclidean space. 

In the next section, we will introduce our assumptions on the functions fi, 
as well as recall Carlier and Ekeland's formulation of the optimal matching 
problem and key theorems from [9] and |29j which we will need for the proof of 
our main result. In the third section, we prove that the solution v to Carlier and 
Ekeland's problem (IMAM I) is absolutely continuous with respect to Lebesgue 
measure. In the fourth section, we prove our main result, while the fifth section 
is reserved for examples and applications. 

2 Preliminaries 

In this section, we present assumptions we will make on the functions /i, and 
recall some key ingredients which will used in the proof. 
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2.1 Definitions and assumptions on the fi. 

Assuming S C'^{Mi), we define D'^,^fi(xi, z) to be the nxn matrix of mixed, 
second order partial derivatives, 

2 f^^c 

Other matrices of second order derivatives will be denoted analogously. We will 
say that fi is non- degenerate if D^,^fi{xi, z) is invertible for all (cc^, z) € Mi x Z. 

We will denote by Dx^fiixi, z) (respectively D^fiixi, z)) the differential of 
fi with respect to Xi (respectively z). We say fi is (a;,, 2;)-twisted (respectively 
(z, a;i)-twisted) if for all Xi £ Mi (respectively z £ Z), the mapping 

z Dx.fi{xi,z) 

(respectively Xi i-> D:,fi{xi,z)) 

is injective. We will always assume the following hypotheses: 

HI For all i, fi is and non-degenerate. 

H2 For each (xi, a;2, Xm) the supremum in ([3|) is attained by a unique z = 

z{xi,x2, e z. 

H3 B(a::i,a;2, ...Xm) := X^Ili -D^z/»(a:^i,^(a;i, ^2, is non-singular. 
H4 /i is (xi, 2;)-twisted. 

Note that condition H2 defines a mapping z : Afi x M2... x M„i — > Z, where 
'z{xi,X2, is the unique maximizer in ([3]). 

In addition, we will sometimes assume 

H5 For all i, fi is (z, a;i)-twisted. 

It was shown in [29], Proposition 3.3.1, that conditions H1-H3 above imply 
that the surplus b is and the following formulae hold: 

Dx.b{xi,X2,-.;Xm) = Dx.fi{xi,z{xi,X2,—;Xm)) (4) 
Dx^z{xi,X2, ■■■,Xm) = -B^^ ■ {xi,X2, ...,X„i)Dl^J^{xi,z{xi,X2, .■■,Xm)) (5) 

Dl^xMixi,X2,...,x^) = (6) 
= -Dl.^fi{xi,z{xi,X2, ....,Xm)) ■ B^^{xi,X2, —,x.m) ' Dl^Jj{xj,z{xi,X2, for i ^ j. 

2.2 Carlier and Ekeland's formulation of the matching 
problem 

Carlier and Ekeland developed an alternative formulation of the optimal match- 
ing problem which will be crucial in our argument [5] ; we briefly review it here. 
Given a Borel probability measure v an define 

TfXv,Hi) svip j fi{xi,z)dn. (7) 

JMixZ 
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where the supremum is over ah measures tt on Mi x Z whose marginals are v 
and fj,i. The muhi-agent matching problem studied by Carlier and Ekeland is 
to find the v which maximizes 

m 

(MAM) 

i=l 

We now record two results of Carlier and Ekeland which we will use in a crucial 
way. First is the existence and uniqueness of a maximizer j/ (Proposition 4 in 

my- 

Theorem 2.2.1. Assume conditions HI and H4 and that /ii is absolutely con- 
tinuous with respect to Lebesgue measure. Then there exists a unique maximizer 
V for (|MAMp . 

The second result of Carlier and Ekeland which we will need encodes a 
natural relationship between the the problems (jMAMp and (IMKp (Proposition 
3 in i9j): 

Theorem 2.2.2. Assume conditions HI, H2 and H4- Then, for the surplus 
b defined by ^ , the following hold. 

l.If^ solves jMQ, then ~ V maximizes (|MAM[) . 

2. If v maximizes (jMAM[) . then there exists a solution 7 to (jMKp such that 
2.3 Rectifiability of the optimizer 

We will also need the following result, which is a special case of Theorem 2.3 in 
PU] : see also Proposition 3.3.1 in [29 for more details. 

Theorem 2.3.1. Assume H1-H3. Let 7 be any minimizer for (jMKp and 

assume conditions H1-H3. Then the support of j is n-rectifiable. In addition, 
it is spacelike for the symmetric bilinear form 

J2[Dl.^b + Dl,^b]. 
That is, for any tangent vector v = {vi,V2, ...,Vm) to spt{'-f), we have: 

i>2 

In view of the preceding result, we will let t £ T , where 

T := spt{-i) C Ml X M2 X ... X Af™, (8) 

and 

Xl{t) G Mi,X2{t) € M2,....,Xm{t) e Mm (9) 
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will denote the coordinate maps. Note that, as xi : T — ?> Mi is certainly Lips- 
chitz, if /ii is absolutely continuous in local coordinates, 7 must be absolutely 
continuous with respect to n - dimensional HausdorfF measure iJ" on T. Fur- 
thermore,, the maps 1 1— >■ Xi{t) are all difFerentiable wherever the rectifiable set T 
is smooth enough to differentiate, which is H"" almost everywhere and therefore 
7 almost everywhere. Together with the spacelike condition in the preceding 
theorem, this implies that 

Y,{DtXlit)f ■ Dl^^fi{xi{t),X2(t), ...,Xra{t)) ■ DtX,(t) > 
i>2 

7 almost everywhere. 

3 Absolute continuity of u. 

In this section, we prove that the solution i' to (|MAMp is absolutely continuous 
with respect to local coordinates. This can be seen as a structural result on i/, 
interesting in its own right, but it is also a key part of the machinery needed to 
prove uniqueness and structural results for (|MKp . 

Lemma 3.1. Assume conditions HI-H4 and that fj,i is absolutely continuous 
with respect to Lebesgue measure. Then, for any solution 7 to (jMK[) . 

7[{i: |det(Aa^i(i))| =0}]=0, 

where Xi{t) is defined by 

Proof. Set A ^ {t : \det{DtXi{t))\ = 0}. We first show that the image of A 
under xi, 

xi{A)^{xi{t) : \DtXi{t)\ =0} 

has Lebesgue measure zero; this is a simple application of the Lipschitz area 
formula. Indeed, the coarea formula tells us that 

I \detDtXi{t)\dH^{t)= [ #[Anx^'{y}]dH^iy)= [ #[A n x^' {y}]dH^ (y) . 

Ja JMi Jxi(A) 

Here, i?" denotes n-dimensional Hausdorff measure on both T and Mi C 
K" (of course, in the latter case, this coincides with Lebesgue measure). As 
\detDtXi{t)\ = on A by definition, the left hand side above is equal to zero. 
On the other hand, #[A n a;^^{y}] > 1 for all y £ xi{A), and so 

/ 4[A^xi'{y}]dH^{y)> I idH^\y)- 

Jxx[A) Jxi(A) 

The right hand side above is simply the Lebesgue measure |a:i(A)| of xi{A), 
and so combining the above two equations yields |xi(j4)| = 0. 



7 



By absolute continuity, iJ,i{xi{A)) = 0. Now, note that 
A^{t: det\DtXi{t)\ = 0} C x^^{xi{A)) = x^^{{xi{t) : det\DtXi{t)\ = 0}), 
and therefore 

7(A) < ^{x^'ix,{A))) 
= 0. 

□ 

Lemma 3.2. Assume conditions H1-H5. Then the mapping T 3 t 'z{xi{t), X2{t), Xm{t)) G 
Z has an invertible derivative almost everywhere, with respect to n-dimensional 
Hausdorff measure on T . 

Proof. Wherever is differentiable (which is almost everywhere, 

by Theorem 1 2. 3. II) we have, by the chain rule and the derivative formulae ([5]): 

m 

Dtz{xi{t),....Xra{t)) ^ '^D^^z{xi{t),...,Xm{t)) ■ DtXi{t) 

i=l 

m 

= -'^[B{xi,X2,—,Xjn)]~^ ■ Dl^Ji{xi,z{xi.X2,—,Xm)) ' DfX^, 
i=l 

where we have suppressed the argument t in the last line. Multiplying by 

[Dtxiit)]'^ ■ [Dl^Ji{xi{t),z{xi{t),X2{t),...,x^{t)))], 
and suppressing the arguments of all the functions, we have 

m 

{Dtx^f ■ Dl^J, ■ D{z = - ^(Axi)^ • Dl^J^ ■ ■ Dlj, ■ D,x. 

m 

= -{Dtxif ■ Dl^J, ■ ■ Dl^Ji ■ Dtxi -Y^Dtxl ■ dI^J^ ■ B^' ■ Dj^J, ■ DtXi 

i=2 
m 

= -{Dtx^f ■ Dl^J, ■ B'^ ■ Dl^J, ■ Axi + J2 Dt^i ■ Dl^.fi ' 

1=2 

(10) 

Note that we have used formula 

Now, by the maximality of z i— >■ X]i=i fi{^ii ^) 'z{xi,X2, Xm), B{xi, X2, Xm) ■= 
J2i=i ^zzfii^ij ^) is negative semi-definite; by assumption H3 it must therefore 
be negative definite. Therefore, at every point where \DtXi{t)\ ^ 0, 

-(Axi)^ • Dl^J, ■ B-' ■ Dl,J, ■ Dtxi 
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is positive definite. As the second term in ([TU| is positive semi-definite by 
Theorem 12.3. 1[ the left hand side 

Dtxl ■ Dl^Ji ■ Dtz 

must be positive definite and hence invertible. Therefore, Diz must be invert- 
ible, as desired. □ 

As we recalled in Theorem 12.2. 1[ Carlier and Ekeland proved that there 
exists a unique optimizer 1/ to (IMAMp under the conditions HI and H4 [9]. We 
prove below the main result of this section; that under conditions H1-H5, this 
optimizer is absolutely continuous in local coordinates. 

Theorem 3.3. Assume fii is absolutely continuous with respect to Lehesgue 
measure. Under conditions H1-H5, the maximizer v to (jMAMp is absolutely 
continuous with respect to Lebesgue measure. 



Proof. By Theorem 12. 2. 2| we have = z#7. Now, Theorem 12.3.11 implies that 

7 =: h{t)dt 

is concentrated on the n-rectifiable set T, and is absolutely continuous with re- 
spect to n-dimensional Hausdorff measure (see the discussion in subsection 2.3); 
we denote its density by h. Furthermore, the map T 9 t z{xi{t),X2{t), ...,Xm{t)) & 
Z is differentiable with an invertible derivative almost everywhere, by Lemma 

Now, choose a set A C Z, negligible in local coordinates. We need to show 
i^iA) = 0, or, equivalently j{z^^{A)) — 0. Let 5 : T — >• R be the characteristic 
function of the and set 

fit) ^ g{t)h{t)/\det{z{xi{t),X2{t), ...,x^{t)))\. 

Then the co-area formula applied to / yields: 

i^iA) = ^iz-\A)) 



git)h{t)dH^{t) 

/ (t ) I detDtz{xi {t),X2{t),...,x,n{t))\dH"{t) 



The first term above is zero because \A\ = 0. We now show that the second 
term is zero as well. 

Now, note that ii z ^ A, then z~^{z} does not intersect 'z~^{A). Therefore, 
for all t £ z^^{z}, g{t) — and therefore, for almost all t S z^^{z}, f{t) = 0. 
It follows that f{t) — for all t G z~^{z}, for z E Z \ A, and therefore the 
integrand in the second term above is zero. This completes the proof. □ 
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4 Monge solutions and uniqueness in the multi- 
marginal problem 

In this section we prove our main result; namely, that the optimizer 7 for (jMKp 
is unique and induces a Monge solution. 

Theorem 4.1. Assume /ii is absolutely continuous with respect to Lehesgue 
measure. Under conditions H1-H5, there exist Borel maps Fi : Z ~> Xi such 
that any optimizer 7 for (|MK[) is equal to the pushforward (Fi, F2, Fm)^i^, 
where v is the unique optimizer for (jMAMp . 

Proof. By Theorem 13.31 v is absolutely continuous with respect to Lebesgue 
measure and by H5, /,; is (z, Xi)-twisted. Note also that as the measure 7 on 
the compact set Mi x M2 x ... x M,„ has compact support and z is continuous, 
u = z^"f has compact support as well. It is a well known result in optimal 
transport theory that these two conditions imply that each optimal coupling 7^ 
in ([7]) is unique and is concentrated on the graph {z, Fi{z)} of a unique mapping 



we will show 7 = 7. 

Carlier and Ekeland proved that for any optimal 7, we have z^7 = v (The- 
orem [2221). They also proved that the measure 

{Tri,z{xi,X2,...,Xm))#J 

on Mi X Z is optimal in where Wi : Mi x M2 x ... x Mm -> Mi is the canonical 
projection. It follows from the remarks above that 

(7ri,z(a;i,a;2,...,x,„))#7 = 7i = {Fi,Id)^v. 

Note that this implies that, 7 almost everywhere, we have 

Xi = TTi{xi,X2, ...,Xjn) = F^{z{xi,X2, ...,Xm))- 

Therefore, 7 almost everywhere we have 

{X1,X2, ...■,Xm) = (^Fi{z{xi,X2, . . . , X„i )) , F2 (z (xi , 2:2 , Xm)) , Fm{z{xi, X2, ...jXm)) 

:= F{z{xi,X2, ...,x„i)), 

where we have defined F : Z ^ Mi x M2 x ... x M„ as F := {Fi,F2, ...,Fm)- 
Now, given any continuous function H : Mi x M2 x ... x Af,„ — > M, this implies 



J Ml X Af2 X...X Af„, 

As e C"(Afi X M2 X ... X Mm) was arbitrary, this implies 7 = 7, as desired. □ 



F^:Z-^ Mi. Set 



j = {Fi,F2,F3,...,Fm)#iy; 
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Corollary 4.2. Under the conditions in Theorem \4-l\ the optimal measure 7 
is unique. 

Proof. The result follows immediately from the uniqueness of v and the maps 
Fi in the preceding theorem. The former is proven by Carlier and Ekeland 
(Theorem 12.2.11 here) . while the latter follows easily from the standard theory 
of optimal transportation. □ 

Finally, we show that the (unique) optimal measure 7 induces a Monge 
solution. 

Corollary 4.3. Under the conditions in Theorem \4.1\ the optimal measure 7 
has the form (Jd, G2, G3, ...Gm)#Mi7 for unique measurable maps Gi : Mi 
Mi. 

Proof. Setting d — FiO F^^ , we easily obtain {Id, 6*2, G3, ...Gm)^/xi =7. □ 

Remark 4.4. (economic interpretation) It seems appropriate to briefly 
discuss the economic interpretation of Theorem \4.1\ in terms of the matching 
models in 111] and I'ri particular, those two papers discuss slightly different 
notions of purity. 

Carlier and Ekeland defined an equilibrium to be pure when workers of the 
same type necessarily signed contracts of the same type. In mathematical terms, 
this means that for the optimal v in (|MAM|) , each optimal coupling 7^ between 
fii and V is concentrated on the graph of a function Fi : Mi Z . 

Chiappori, McCann and Nesheim, on the other hand, defined an equilibrium 
in the m = 2 case to be pure when buyers of the same type almost always 
conducted business with sellers of the same type. This means that the optimal 
measure 7 in (jMKp is concentrated on the graph of a function G2 '■ Mi M2 • 

The purity notion of Chiappori, McCann and Nesheim extends naturally to 
the m > 3 setting. We will say that an equilibrium is pure, in the sense of 
Chiappori, McCann and Nesheim, if the optimal 7 in (jMK[) is concentrated on 
the graph of a function over xi. Economically, this means that workers in Mi 
of the same type almost always conduct business with a team of the same type. 

Purity of the equilibrium in the sense of Carlier and Ekeland is proven in f^; 
Corollary \4-l\ can be interpreted as purity in the sense of Chiappori, McCann 
and Nesheim. 

5 Examples and applications 

In this section we present several examples and applications of our main the- 
orem. The examples in the first subsection demonstrates that Theorem 14.11 is 
a generalization of the seminal paper of Gangbo and Swi§ch |17| : in fact, we 
show that our work here implies a partial generalization of their result to mani- 
folds. The example in subsection 5.2 shows that Corollary 14. 31 also encompasses 
a result of Heinich [18) . In subsection 5.3, we present an example of a surplus 
function which is not of the form ([3]), but which satisfies different conditions 
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required for uniqueness and Monge solutions in [5S]; this verifies that Corollar- 
ies 14.21 and 14.31 do not generalize the main theorem in [28] . On the other hand, 
we exhibit an example in subsection 5.4 which is of the form ([3]), but does not 
satisfy the conditions in [35]. Therefore, the main theorem in does not 
generalize Corollaries 14.21 and 14.31 either. 

5.1 Generalization of the two marginal quadratic cost on 
manifolds 

In this subsection, we consider the case where each Mi C M is a bounded, open 
subset of some common Riemannian manifold M , Z = M and the —fi are all 
the Riemannian distance squared. In this case, we have: 

m m 
b{xi,X2,-.,Xm) = SUpy2-(f{Xi,z) = sup y2-(f{Xi,z) (11) 

Note that maximizing (jMK[) is equivalent to minimizing 

/ C{X1,X2, ...,Xm)dj{xi,X2, ■■■,Xm) 

J A/ixM2X...xMm 

over the same set of measures for the cost function 

m 

c{xi,X2,—,Xm) ^ -b{xi,X2,.-,Xm) ^ inf rf^ (a;^ , z) (12) 

z^M ^ — ^ 

1=1 

Recall that M is called a Hadamard manifold if it is complete, simply connected 
and has non-positive sectional curvature. Our main result easily implies the 
following: 

Corollary 5.1.1. Let M be a Hadamard manifold. Let ^i, /i2, .../!,„ be com- 
pactly supported Borel probability measures on M and assume /ii is absolutely 
continuous with respect to local coordinates. Then the solution 7 to (|MKp with 
surplus b given by (jlll) is unique and is concentrated on a graph over xi . 

Proof. By the Cartan-Hadamard theorem, M is diffeomorphic to M" and has 
no cut locus. Conditions H1,H4 and H5 are then easily verified. Furthermore, 
it is well known that each mapping 

z !—> d'^{xi, z) 

is uniformly convex along geodesies. Therefore, 

m 

y^^d^{xi,z) 

i=l 

is uniformly convex along geodesies, which implies conditions H2 and H3. The 
result now follows directly from Theorem 14.11 and Corollaries 14.21 and 14.31 □ 
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Remark 5.1.2. Gangho and Switch proved uniqueness and existence of Monge 
solutions for the quadratic surplus ^ on Euclidean space; a natural open ques- 
tion is whether this result can he generalized to Riemannian manifolds. Perhaps 
the most obvious generalization of this result would consist of uniqueness and 
Monge solution results for the cost function d'^ixi,Xj)'^ . However, the pre- 
ceding Corollary can be interpreted as a different type of generalization. Note 
that 

m ^ 

sup E-|2^«-^l^ = -El^' ~ -E^jI^ 

i 

z — ^ rn ^ — ^ 



i=l i j 



m 



E(2---E 

z — ^ rn ^ — ^ 



m 

EE 



Xi 
J 



— TTl 





1 


lE-.H 

3 




^ 1 
m 


lE-.i' 

3 




^ 1 
m 


EE^j 

3 k 


Ei^* 

i 













^ Xj \Xi\ ) ___/,/ ,Xj ■ Xk 
3 



This establishes two facts: 



1. The quadratic surplus ([T]) is of the form for fi{xi,z) = —\xi — zp, 
which, as can easily be checked, satisfy H1-H5. Corollaries \4.S\ and \4.3\ 
therefore generalize the theorem of Cangbo and Swi§ch. 

2. In particular, the preceding Corollary shows that our result implies a gen- 
eralization of Cangbo and Swi§ch's theorem to Hadamard manifolds, 

Remark 5.1.3. Brenier's celebrated polar factorization theorem amounts to 
uniqueness and existence of Monge solutions for the optimal transport problem 
on Ml = M2 = R" with surplus — \x — J^. McCann's generalization of this 
result resolves the same questions for the quadratic distance squared —d'^{xi,X2) 
on a Riemannian manifold J23f . 
Note that 



2 

and that on a Riemannian manifold 

^^^^=sM~d{x,zf~d{y,zf). 

Therefore, the preceding Corollary can be seen as a multi-marginal version of 
McCann's theorem (albeit only for Hadamard manifolds). Furthermore, by the 
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previous remark, Corollary \5.1.1\ generalizes the Gangho and Switch result in 
the same spirit that McCann's theorem generalizes Brenier's theorem. 

To summarize the last two remarks: the surplus' of Brenier, McCann and 
Gangbo and Switch are all of the form Q with quadratic functions fi. For 
Hadamard manifolds, then, Corollary 15 .1.11 can therefore be thought of as gen- 
eralizing McCann's surplus to several marginals, or the Gangbo-Swi§ch sur- 
plus to Riemannian manifolds. One should note that, unlike in [23] , we deal 
here exclusively with a class of manifolds which are diffeomorphic to M" and 
therefore avoid entirely complications arising from the cut locus, including non- 
smoothness of the fi and non- uniqueness of the minimizing 'z{xi,X2, Xm)- On 
the other hand, to this point very little is known about multi-marginal problems 
on any non-flat manifolds with surplus functions derived from the distance, and 
so I feel Corollarv lS.l.ll is an interesting contribution in this direction. 

Finally, recall that the barycenter (with equal weights) of several points 
Pi,P2, ■■■,Pm in a metric space is the minimizer of g i— >■ X^I^Li '^^(Pij 9)j assuming 
this minimum exists uniquely. Recall that Wasserstein space over a Riemannian 
manifold M, the space of Borel probability measures on M endowed with the 
Wasserstein metric 



H^^(Mi,z.):=./inf / ^-^dn 

where the infimum is over all measure tt on whose marginals are /i and i', 
is a metric space. 

For the surplus function (fTTj). the optimizer v in (jMAMp in this case coin- 
cides with the barycenter of the measures /ii,/i2, ■■■fJ-m- Barycenters in Wasser- 
stein space over K" were studied by Carlier and Agueh, who proved existence, 
uniqueness and regularity results. The present author proved existence and 
uniqueness of barycenters in Wasserstein space over Riemannian manifolds in 
[25] : the following result follows immediately from Theorem l3.3l and can be seen 
as a generalization of the regularity result of Agueh and Carlier. 

Corollary 5.1.4. Let fii, fi2, ■■■fJ-m be compactly supported Borel probability 
measures on a Hadamard manifold M and assume fii is absolutely continuous 
with respect to local coordinates. Then the barycenter v is absolutely continuous 
with respect to local coordinates. 



5.2 Convex functions of the sum 

Let b{xi, X2, Xm) = ^(X]"=i where h : M" — !■ K is smooth and uniformly 
convex. Letting h* denote the Legendre transform, we have that 

m m 

Xi) = sup y Xi ■ Z — h*{z) 

Letting fi{xi,z) — xi ■ z — h*{z) and fi{xi,z) = Xi ■ z for i > 2, we have 
that b is of the form ([3]). It is straightforward to verify the conditions H1-H5, 
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and so our main theorem generalizes the result of Heinich [TS] , at least under 
smoothness and uniform convexity conditions (the result of Heinich required 
only strict convexity of h) . 

5.3 A non hedonic pricing surplus 

In the last two subsections, we have seen that our main result generalizes the 
results of Gangbo and Swi§ch and Heinich. It is natural to wonder whether it 
also encompasses the result in [5S]. Below, we answer this question negatively, 
by exhibiting a surplus function satisfying the conditions in |28| which is not of 
the form 

Let m — S and set 

b{xi,X2,X3) = xi • a;2 + xi • ^3 + a;2 • Axs (13) 

where A is an n x n matrix which is positive definite but not symmetric. It was 
verified in [28] (Example 4.2) that this surplus function satisfies the conditions 
therein. We now show that this surplus function cannot be of the form ([3]) . 

Proposition 5.3.1. The surplus b defined by ((T3)l is not of the form ([3]) for 
any functions /i, /2 and /a satisfying H1-H6. 

Proof. The proof is by contradiction; assume there exist functions /i , /2 and /a 
satisfying H1-H6 so that b — sup^g^ Y^i=i fii^i^z). We compute the following 
product of matrices of mixed partials: 

On one hand, by (fT3|) . this is 

Dl...b-[Dl,,br'-Dl,^_b^A. 

On the other hand, using (|6]) and suppressing the arguments of the functions, 
we have 

Dl.sb-[Dl.sbr'-Dl,^b 

= -Dlj2 ■ B-^ ■ Dl^Js ■ [Dl^Js]-' ■ B ■ [DlJi]'' ■ DlJ, ■ B'^ ■ D^J^ 
= -Dlzf2 ■ B-' ■ Dl^j2 

We therefore have 

A^-Dlj2-B-'-Dl^j2 (14) 

As B-^ is symmetric and D^.^ [Dl^j2f, ~Dl^j2-B-^-Dl^j2 is symmet- 
ric. On the other hand A is not symmetric, and so ([T4| yields a contradiction. 
This completes the proof. □ 
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5.4 Suprema of convex functions 

We now consider the case where each fi is a concave function of the sum. 

Corollary 5.4.1. Assume for all i that fi{xi,z) = hi{xi + z) for a smooth, 
uniformly concave function hi such that \yai]^y\^oo hi{y) — — oo and that the 
surplus is given by ([3]). // /ii is absolutely continuous with respect to Lebesgue 
measure, then the optimal measure 7 for (|MKp induces a Monge solution and 
is unique. 

Proof. We need only to verify conditions H1-H5, in order to apply Corollaries 
I4.2l and [4.3l HI, H3, H4 and H5 follow from the smoothness and uniform con- 
vexity of hi, while the growth condition ensures the attainment of a maximum 
(condition H2). Therefore, our main theorem applies. □ 

Recall that in [28], a strong second order condition on the surplus, called 
condition (III) in Theorem 3.1 there, is required to prove Monge solution and 
uniqueness results. Next, we show that, for a surplus of the form (jSj, for well 
chosen fi{xi, z) = hi{xi + z) (satisfying the hypotheses in the preceding Corol- 
lary) , the condition (III) fails, and so our result here is not strictly weaker than 
the theorem in |28j . 

For simplicity, choose to = 3 and take each to be a large open ball, 
centered at the origin. In [35], condition (III) was calculated explicitly for 
surpluses of form Q in Proposition 4.2.2; amounts to the following: 

For all choices Xi,xi € Mi, X2 G M2 and X3,X3 € M3, the n x n matrix 



{Dx2zf2)B ^^^2/2) {x2,z{xi,X2,Xs.)^ + Dl,^^^f2{x2,z{xi,X2,X:S) - 01^X2^(^2, z{xi, X2, X3) 



is positive definite. 

We now show this condition fails for the choice fi{xi,z) = — y^T+|a!7+2p] , i 
1, 2, 3. This form of preference function was introduced in the optimal transport 
literature by Brenicr [5^ . It also served as an early example of a cost function 
satisfying (A3w), the crucial condition introduced by Ma, Trudinger and Wang 
governing the regularity of optimal maps in two marginal problems [22) . 

It is easy to check that this function satisfies the conditions in the preceding 
Corollary, however, and so our surplus ([3]) with these fi is an example of a 
surplus to which our main result here applies but the result in [28| does not. 

Proposition 5.4.2. Condition (III) above fails for b of the form when 
TO = 3 and each fi{xi, z) = — \/[l -I- \xi + zp]. 



Proof. Take 



Xl,X2,X'i = U 2:3 =Xi 



for some p € M" with \p\ = |. As DJ,{x,,z) = Dh,{xi + z) ^ 



we then have 



[l+\xi+z\^]^' 
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= ^D^fi{Xi,z{Xi,X2,X3)) 



i=l 

3 



^m(0 + z(0,0,0)) 



i=l 



3^(0,0,0) 



[1 + |^(0,0,0)|2]^ 



which implies 



z{Xi,X2,X3) = z(0,0,0) = 0. 
In a similar manner, we have 

= Dzfi{xi,z{Xi,X2,X3)) + Dzf2{x2,z{xi,X2,X3)) + L'z /a (X3 , Z (f i , ^2 , fa ) ) 

= Dh{p + z{p, 0, p)) + Dh{0 + zip. 0.p)) + Dh{p + z{p, 0, p)) 

^ p+-{pA),p) np-i),p) 

[l + \p + z{p,0,p)\^]i [l + \z{p,0,p)\^]i 
One can easily verify (by, for example, direct substitution) that this yields 

4 

z{p,0,p) = --p. 
In particular, as |p| = |, 0,p)| = 2. Now note that 



J ( y-y'^ 



D'h{y) = j^^^{^rh2-I)- 



It is then straightforward to compute 

B{xi,X2,X3) = -B(0,0,0) = 3Z)''(0 + 0) = -3/ 
Dl^zh{x2,z{xuX2,X3) = D^h{0 + 0) = -I 

Therefore 

\Dl,j2)B-'iDl^j2)\ (^X2,z{xi,X2,X3)) = [(1?^, Jz)^^^ (i^L. /s)] (o, o) 

Similarly, we have 
^S2X2/2(2^2,^(a;i,a;2,a;3)) = -/, 

D.^.J.(x2,zix,,X2,X3)) = [i + |^(^,o,p)P]i( l + mP,0,pW 



17 



We therefore have: 



rp _ \t_t \ ^ Z(p,0,p) ■ Z(p,0,p)^ 

3 [1 + |-(p,0.p)|2]2 1 + |z(p,0,p)|2 

2 . 1 
< -7:1 + 



3 [l + |^(p,0,p)|2]i 
[ , - V 



As < |, this is negative definite; in particular, it is not positive definite, so 
condition (III) fails. □ 
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